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A  On*-Dlm*n*lonal  Vertical  Diffusion  Paromatar 

for  Extramaly  Inhomoganaout, 
Layarad  Turbulanca  in  Stratlfiad  Fluid* 


l,  INTRODUCTION 

The  mntn  motivation  for  this  Investigation  la  the  preaent  concern  about  strato¬ 
spheric  pollution.  The  stratosphere  Is  characterized  by  a  high  degree  of  stability, 
and  turbulence  In  this  region  occurs  In  the  form  of  thin,  random,  horizontal, 
pancake -shaped  layers.  *  The  shear,  or  Kelvln-Helmholtz,  Instability  Is  presumed 
to  be  'he  mechanism  for  the  formation  of  these  layers.  The  nature  of  this  type  of 
turbulence  ta  discussed  In  Rosenberg  and  Dewan,  1  Dewan  <  1079a),  3  and  Dewan 
(1079b),3  There  are  two  purposes  for  the  present  report.  The  first  Is  to  overcome 
some  of  the  Inltlrl  criticisms  we  received  (Rosenberg  and  Dewan1)  due  to  the  fact 
that  we  used  the  term  "eddy  dlffuslvlty,"  As  will  be  seen,  we  will  be  considering 
an  entirely  new  parameter  which  Is,  In  a  sense,  completely  different  from  the 
usual  eddy  diffusion,  This  new  parameter,  like  eddy  diffusion,  Is  phenomenological 
In  nature;  however,  this  Is  where  the  similarity  ends.  To  explain  the  new  param¬ 
eter,  this  report  shall  make  use  of  an  analogy  between  the  random  layer  mixing 

(Received  for  publication  13  June  1980) 
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process  mentioned  above  and  the  finite  difference  formalism  used  to  numerically 

simulate  the  diffusion  equation,  That  such  an  analogy  exists  may  surprise  some 
readers!  nevertheless,  It  represents  a  good  way  to  convey  the  physios  of  the  proc¬ 
ess,  This  report's  second  purpose  Is  to  demonstrate  that  our  new  parameter  la 
Indeed  a  self-consistent,  "diffusion  like"  parameter,  Thto  will  be  achieved  by 
means  of  a  series  of  digital  computer  simulations, 

The  usual  eddy  diffusion  parameter  (eddy  diffusion  coefficient)  for  turbulence 
transport  assumes  that,  at  least  as  a  practical  approximation,  the  turbulence  Is 
essentially  homogeneous  in  some  sense,  Paaqulll*  defines  eddy  dtffustvity,  K,  as 

K*  vi  (1) 

where  v  and  i  are  the  velocity  and  size  scale  lengths,  respectively,  which  are 
appropriate  for  the  turbulence  in  question.  Another  form  to  be  found  In  the  litera¬ 
ture  (Lilly  et  al,  5  Panofsky  and  Heck8)  la 


K 


(2) 


where  T  U  the  average  eddy  dissipation  rate  and  Nfi  la  the  buoyancy  frequency. 

1,1  Assumption*  and  Definition  of  “Bulk  Diffuiivity" 

We  shall  assume  that  (a)  the  turbulent  layers  occur  at  random  heights  and  at 
random  ttmss  (this  assumption  ts  usually  made  in  the  context  or  stable  riulds, 
see  for  example,  Panofsky  and  Heck,  8  Lllley  et  al,  5  Woods,  7  and  Woods  and 
Wiley8),  (b)  there  is  no  vertical  transport  outside  of  the  turbulent  layers  (In  other 
words,  between  layers,  motion  is  regarded  as  laminar),  and  (c)  within  the  mixing 
layers,  the  mixing  Is  presumed  to  be  total,  Regarding  this  last  assumption,  It 
remains  to  bo  proven  correcti  on  the  other  hand,  certain  experimental  observations 


4.  Pasqulll,  F.  (1C74)  AtmOpherlc  Diffusion.  2nd  Ed.,  John  Wiley  &  Sons. 

8,  Lilly,  D.  K. ,  Waco,  D.  E.,  and  Adelfang,  S,  I.  (1875)  Stratospheric  mixing 
estimated  from  high  altitude  turbulence  measurements  by  using  energy 
budget  techniques,  The  Natural  Stratosphere  of  1874.  CIA P  Monograph  1, 
Final  Report,  DOT-TST-75-51,  pp,  8-81  to  Q-60, 

0.  Panofalw,  H.A,  and  Heck,  W.  (1975)  Stratospheric  mixing  estimates  from 

SHU1 t  C,AP  M°no- 

7.  Woods,  J.D,  (1968)  Wave-Induced  shear  instability  in  the  summer  thermo- 

cline,  J,  Fluid  Much,  32t7Bl-B00, 

8.  Woods,  J.D,  and  Wiley,  R,  L,  (1072)  Billow  turbulence  and  ocean  micro¬ 

structure,  Deep  Sea  Research  and  Oceanic  Abst.  10rB7-121. 
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have  been  made  which  are  consistent  with  it.  These  are  the  temperature-profile 

g 

studies  of  Mantis  and  Peppln,  who  note  that  balloon-borne  sensors  indicate  that 
in  the  lower  stratosphere  and  upper  troposphere  there  are  regions  of  nearly 
adiabatic  lapse  rates,  Analogously,  in  the  upper  ocean  where  the  stability  situa¬ 
tion  is  very  much  like  the  stratpsphere,  Woods  and  Wiley8  have  reported  step- 
like  structures  in  the  vertical  temperature  profiles.  Such  profiles  are  consistent 
with  the  hypothesis  that  total  mixing  can  result  from  the  turbulent  layers,  Finally, 
one  more  assumption  is  needed)  (d)  there  is  enough  horizontal  homogeneity  so  that 
a  one-dimensional  vertical  transport  model  will  be  adequate, 

In  the  next  section,  the  mechanism  for  vertical  transport  by  means  of  random 
mixing  layers  will  be  explained.  At  this  stage,  however,  the  assumption  that  total 
mixing  takes  place  can  be  used  to  suggest  (by  means  of  a  dimensional  argument) 
that  any  effective  "bulk  diffusion  coefficient, "  which  parameterizes  layered  turbu¬ 
lence,  could  not  be  appropriately  described  by  Eqa.  (1)  or  (2),  The  argument  la 
based  on  the  fact  that  the  width  of  the  layer  and  the  time  intervals  between  layer 
formations  would  determine  the  transport.  These  parameters  do  not  explicitly 
appear  in  Eqs.  (1)  or  (2),  and  unless  there  were  a  proven  connection,  say,  between 
Ng  and  the  time  Interval  between  new  random  layer  profiles,  it  Is  not  at  all  clear 
that  Eqs,  (1)  or  (2)  could  be  useful  in  the  present  context,  on  the  other  hand,  the 
parameter  to  be  developed  here  will  explicitly  depend  on  layer  thickness,  A,  and 
the  relevant  time.  Interval,  At. 

Before  proceeding  further,  it  la  essential  to  precisely  define  Kg,  the  bulk 
vertical  diffuslvlty  parameter.  Let  F  designate  the  average  flux  of  material 
through  a  given  slab  of  stable  fluid.  This  slab  is  assumed  to  be  very  much  larger 
than  the  largest  turbulent  layers,  Let  4  designate  a  scalar  quantity,  such  as  poten¬ 
tial  temperature  or  mixing  ratio  of  some  neutrally  buoyant  substance  (for  example, 
pollutant),  Then 

Kg  .-JL.  (3) 

U  8?/ 07. 

where  the  overbar  designates  an  average,  This  definition  is  based  on  Fourier's 

2  fl 

heat  transport  equation  (cf.  Dewan  as  well  as  Panofsky  and  Heck  ). 


9.  Mantis,  M.T,  and  Pepin,  T.J.  (1071)  Vertical  temperature  structure  of  the 
free  atmosphere  at  mesoscale,  ,T.  Oeophys,  Res,  30:803 1-8628, 
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2.  THE  MECHANISM  OF  TRANSPORT  BY  MIXING  LAYERS 


Figure  1  depicts  a  slab  of  stable  atmosphere,  and  within  the  slab  is  located  a 
single  turbulent  layer.  As  previously  mentioned,  we  assume  that  no  vertical 
transport  exists  outside  of  such  a  turbulent  layer,  but  let  us  temporarily  add  one 
more  assumptions  (e)  the  layer  remains  fixed  at  one  altitude.  We  ask  the  ques¬ 
tions  "What  is  the  flux,  F,  between  points  Y  and  Z,  located  at  the  top  and  bottom 
of  the  slab,  respectively?"  Obviously,  there  will  be  no  transport  of  material  from 
point  Y  to  Z,  Equation  (3)  therefore  implies  Kg  >  0,  since  there  F  Is  to  be  inter¬ 
preted  as  a  flux  from  poLnt  Y  to  point  Z,  We  have  an  analogy  here  with  the  situa¬ 
tion  in  which  there  la  an  electrical  conductor  located  between  two  regions  of  insula¬ 
tion.  One  cannot  average  the  conductivity  of  the  "sandwich,"  One  is  left  with  a 
nonconductor  (In  the  vertical),  no  matter  how  low  the  resistance  of  the  conductor 
inaide  may  be.  These  considerations  Imply  that  the  intensity  of  turbulence  (once 
total  mixing  is  assumed)  is  of  no  direct  relevance,  and  that,  as  will  be  made  more 
clear  below,  the  temporal  behavior  of  the  layers  is  absolutely  crucial  for  transport 
effects. 

Next  consider  Figure  2,  Here  we  consider  a  temporal  sequence  of  Layer  for¬ 
mation.  Layer  I  is  the  first  to  form.  Total  mixing  takes  place,  and  then  the 
turbulence  decays  In  time,  Due  to  the  gradient  of  4  which  is  indicated,  the  mixing 
causes  material  In  the  upper  half  (A)  of  the  layer  to  move  into  the  lower  half  (B)  of 
the  layer.  Next,  imagine  Layer  II  to  form,  and  notice  that  Layer  II  significantly 
overlaps  the  original  position  taken  previously  by  I.  When  II  has  mixed,  the  profile 
of  4(z)  would  again  be  altered  to  a  single  average  value  within  the  layer.  This  Is 
due  to  the  assumption  of  total  mixing,  This  process  would  result  In  a  further  net 
transfer  or  d  (for  example,  material)  downward.  It  should  be  clear  that  Some 
material  which  was  initially  near  the  top  of  Layer  I  (for  example,  in  region  A)  can 
find  Its  way,  at  the  end  of  this  two-fold  mixing  process,  to  a  point  near  the  bottom 
of  Layer  II.  This  simple  example  demonstrates  the  essential  mechanism  for  ver¬ 
tical  transport  by  means  of  a  random  sequence  of  mixing  layers  In  an  otherwise 
stable  fluid  (all  other  flow  being  laminar  and  horizontal).  Such  a  process  can 
eventually  give  a  net  flux  from  the  top  to  the  bottom  of  the  slab.  Notice  that  there 
are  two  conditions  needed  for  the  mechanism  to  functtom  (a)  a  large  amount  of 
mixing  within  the  layer,  and  (b)  an  overlap  of  layer  positions. 

After  the  random  formation  of  turbulent  layers  In  the  above  fashion,  one  could 

calculate  an  average  flux  F  across  the  slab,  as  well  as  8d/8z.  could  then  be 

1  o  *•* 

calculated  from  Eq.  (3).  In  Rosenberg  and  Dewan  and  Dewan,  It  was  shown 
that  one  can  calculate  an  estimate  for  (designated  here  as  KB.CAI  c)  If  one 
assumes,  among  other  things,  a  steady  state,  This  will  be  discussed  in  the 
Appendix,  The  result  Is 
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Figure  1.  Single  Stationery  Mixing  Layer  Within  »  Slab  of  Atmosphere.  No 

transport  In  this  case 


Figure  2,  Mechanism  of  Vertical  Transport  by  Sequence  of  overlapping  Mixing 
Layers  In  the  Presence  of  a  Gradient 


(4) 


^■CAUC 

where  A  la  layer  thtckneaa,  the  overbar  la  an  average,  and  la  the  average 
time  Interval  between  layer  formattona  overlapping  the  bottom  boundary  of  the 
slab  (or  any  preaet  altitude,  far  that  matter).  Equation  (4)  la  baaed  on  Eq,  (3). 


3.  RANDOM  MIXING  AS  A  STOCHASTIC  STIMULATION  OF  THE 
DIFFUSION  EQUATION 

The  dlffualvlty  given  by  Eq.  (4)  la  new  (aa  of  Rosenberg  and  Dewan1),  and 
before  It  can  be  fully  accepted  aa  valid,  It  muat  be  ahown  to  be  aelf-conalstent,  In 
other  worda,  the  question  remains  as  to  whether  or  not  random  mixing  layer  trans¬ 
port,  aa  estimated  from  Eq,  (4),  la  dlffualve  In  nature  In  the  aenae  that  one  could 
use  the  same  Kg  to  deacrlbe  the  transport  effecta  on  any  Initial  distribution  of 
<Mz).  In  addition,  It  Is  essential  to  establish  Kg  on  a  more  reliable  physical/ 
mathematical  basts  than  has  been  done  heretofore.  The  remainder  of  this  paper 
will  be  devoted  to  these  objectives,  and  the  purpose  of  the  present  section  Is  to 
connect  the  random  mixing  layer  process  to  the  finite  difference  formulation  of  the 
diffusion  equation.  (This  analogy  was  mentioned  previously  In  the  Introduction. ) 

It  la  well  known  that  If  one  wishes  to  stmulate  a  Laplaclan,  by  means  of 

a  finite-difference  numerical  scheme,  one  averagea  a  given  lattice  point  with  all 
of  its  nearest  neighbors,  as:  replaces  the  original  value  at  the  given  lattice  point 
with  tnia  average.  This  process  la  repeated  for  each  time  step.  Now  consider  the 
eflect  of  a  mixing  layer.  In  effect,  the  layer  averages  the  values  of  $(z)  within  Its 
boundaries  and  replaces  these  original  values  with  the  average  (see  Figures  1  and  2), 
This  leads  to  the  Idea  that  random  turbulent  layers  In  a  stable  fluid  could  stochas¬ 
tically  simulate  the  one-dlmenalonal  diffusion  equation. 

In  order  to  make  this  Idea  more  concrete,  we  first  consider  a  specific  pattern 
of  a  mixing  layer  sequence  which  results  In  an  exact  simulation  of  the  finite  differ¬ 
ence  diffusion  equation.  This  pattern  will,  of  course,  be  extremely  regular,  and 
In  this  respect,  will  differ  from  the  random  situation  postulated  for  stratified  tur¬ 
bulence.  Nevertheless,  it  will  serve  to  establish  an  Important  connection  between 
the  two  processes. 

Figure  3  shows,  In  epoch  l,  a  vertical  set  of  layers,  each  In  sequence  immed¬ 
iately  above  the  other,  They  nearly  touch,  but  there  Is  no  overlap.  These  are  to 
be  Imagined  as  mixing  until  all  pairs  of  altitude  points  are  averaged  (we  assume 
each  layer  contains  two  altitude  grid  points).  In  principle,  these  layers  would  not 
have  to  mix  simultaneously.  We  assume,  however,  that  the  mixtng  process  has 


Figure  3.  A  "Mixing  Layer  Simulation"  of  the  Fi¬ 
nite  Difference  Diffusion  Equation 


terminated  at  a  apeclflc  time  which  deilgnatea  the  end  ofEpooh  1.  After  this,  the 
original  layera  are  replaced  by  a  different  set  of  layera  (Time  2  In  Figure  3), 
which  preclaely  and  aymmetrlcally  overlap  the  poaltlona  of  the  original  layera.  In 
other  worda,  the  "cracks"  between  the  original  layer  boundarlea  are  now  In  the 
exact  centers  of  the  Epoch  2  layers,  These  layers  then  mix  and  perform  another 
averaging  of  allltude  pairs. 

The  above  double  sequence  Is  then  repreated  for  each  time  step,  and  so  on. 

We  now  calculate  Its  effects.  Let  z,  as  usual,  designate  the  altitude,  and  let  Az 
be  the  distance  between  the  vertical  points,  Let  t  designate  the  time,  and  let  At 
be  the  Interval  of  time  between  steps.  Each  epoch  described  above  corresponds 
to  only  one-half  of  a  time  step,  At/2.  Thus,  for  $,  at  the  end  of  Epoch  1,  we  have 


In  other  words,  at  z,  at  the  end  of  the  first  epoch,  equals  the  average  of  4  taken 
at  z  and  z  +  Az.  The  mixing  layer  Includes  both  z  and  z  +  Az,  the  point  z  being 
located  a  distance  (l/2MAz)  from  the  bottom  of  the  layer,  while  (z  +  Az)  would  be 
<3/2)(Az)  up  from  the  bottom  (layer  thickness  ■  2  Az),  Next,  Epoch  2  of  the  time 
step  (Epoch  2  In  Figure  3)  results  In 


^(t  +  At,  z) 


(6) 
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In  other  words,  this  time  tin-  uverafie,  or  mixing,  involves  the  grid  point  immed¬ 
iately  below  z  (for  example,  z  -  Az).  We  now  wish  to  relate  4  (t  +  At,  z)  to  the 
original  value*  of  at  t.  Therefore,  insert  Eq.  (9)*  into  Eq,  (0)  to  obtain 


^(t  +  At,  z)  «  7  (2*(t,  »,)  +  $(t,  z  +  Ai)  +  d(t,  z  •  Az)] 

4 

Next,  we  put  Eq.  (7)  into  the  form  of  the  finite  difference  equation  a*  follow*) 
flr*t,  the  finite  time  derivative  will  be  defined  by 


.  Mil  -  ♦<*. 


■nd  the  second  spatial  derivative  (*ee  any  book  on  numerical  analyst*)  will  be 
defined 


AjL\t.*l  ,  J &JLt 

Aa* 


la)  -  24(t,  a)  +  4(t,  a  -  Aa) 


Finally,  we  define  K  as 


K  ,  I  (*«)' 


Rearranglng  Eq.  (7),  and  using  Eqs.  (B),  <B)  and  (10),  we  obtain 


#  ,  k-A-4 

At  (Aa)2 


which  la,  of  oourse,  the  finite  difference  diffusion  equation,  (Note  the  resemblanoe 
between  Eqs.  (10)  and  (4).)  Thus,  the  sequence  in  Figure  3  exactly  simulates 
Eq,  (11).  Actually,  there  are  many  more  sequences  whtch  would  effectively  d(  he 
same  thing,  and  the  uniqueness  of  the  process  In  Figure  3  la  merely  Its  simplicity. 

The  above  discussion  carries  an  Important  message)  namely,  that  In  spite  of 
the  extremely  Inhomogeneous  and  discrete  nature  of  the  process,  a  "diffusion" 
parameter,  called  K  above,,  haa  an  "exact"  connection  with  the  process  of  diffusion, 
Layered  turbulence  Is  extremely  Inhomogeneous  In  the  vertical  direction.  It  Is 
difficult  to  Imagine  a  form  of  turbulence  which  Is  more  Inhomogeneous  In  tjils 

"'Equation  (5)  must  be  modified  so  that  It  can  be  used  to  be  substituted  Into  both 
terms  In  the  numerator  of  Eq,  (9), 


respect.  Nevertheless,  we  find  here  a  way  to  meaningfully  assign  n  dlffustvity  to 
the  process. 

As  was  mentioned,  layered  turbulence  in  actual  fluids  would  differ  from  the 
above  process  by  virtue  of  its  stochastic  nature)  layers  would  presumably  form 
at  random  altitudes,  at  random  times  and  with  random  thickness.  This  leads  to 
the  next  question,  Could  such  a  random  proceaa  effectively  almulate  diffusion  In 
a  manner  analogous  to  the  above?  Also,  would  it  do  so  in  a  consistent  manner? 

In  order  to  show  that  the  anawera  to  the  above  questions  are  affirmative,  we 
shall  make  several  digital  simulations  of  the  stochastic  mixing  process.  The 
effective  values  of  Kg  will  be  determined  directly  from  the  simulated  results  by 
means  of  a  comparison  with  continuous  diffusion  theory.  We  shall  see  that  (a)  the 
behavior  of  4(b)  in  time  la  extremely  similar  to  the  predictions  of  continuous 
theory,  and  (b)  the  Kg  values  Inferred  from  the  behavior  of  the  simulations  are  in 
surprisingly  accurate  quantitative  agreement  with  one  another. 

4,  SIMULATION  PROGRAM  FOR  Kg  »  F7<64/0t) 

In  the  following,  a  digital  program  for  simulating  the  effects  of  random  mixing 
layers  will  be  described.  More  specifically,  a  program  for  ascertaining  the  value 
of  Kg  on  the  basis  of  Eq.  (3)  (see  this  section  title)  will  be  discussed  in  some 
detail,  With  minor  modifications,  the  same  program  also  provides  the  simulations 
to  be  given  In  Section  G.  At  the  first  reading,  some  readers  may  wish  to  skip  some 
of  the  details. 

Figure  4  shows  tho  essential  aspects  of  how  the  simulation  was  done.  A  "alab" 
of  thickness  R  is  depicted  In  an  "environment"  (as  *  0  to  z  =  max).  This  slab  could, 
for  example,  represent  the  entire  stratosphere  or  upper  ocean,  or  it  could  repre¬ 
sent  a  smaller  height  region  of  interest,  It  Is  the  region  for  which  we  seek  a  bulk 
value  of  Ulffuslvlty.  The  reason  for  the  "environment"  will  be  explained  below. 
Inside  the  slab  is  a  "turbulent"  layer  of  random  thickness  A  (at  random  altitude,  as), 
and  within  such  a  layer,  the  values  of  4  would  be  replaced  by  the  average  value 
within  the  layer.  The  Initial  value  of  4  la  shown,  i.nd  It  is  a  straight  line  with  4 
differing  by  unity  In  going  from  the  bottom,  B,  to  the  top,  A,  or  the  slab,  In  other 
words,  the  slope  of  4  across  the  slab  is  (1/R),  as  shown)  the  same  slope  continues 
throughout  the  slab’s  "environment". 

The  height  z  is  chosen  by  means  of  a  random-number  generator.  As  for  the 
thickness,  A,  the  following  procedure  was  used.  A  Becond  random  number  gen¬ 
erator  (which  was  Independent  from  the  one  used  for  z)  was  used  to  pick  numbers 
from  0  to  98,  The  possible  thickness  of  layers  were  3Az,  5A*,  7Az,  DAz  and 
UAz,  that  la,  the  layers  were  to  contain  from  3  to  11  points  (odd  only),  and  In  the 
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UPPER  CURSOR 


calculations  As  -  50  m.  The  use  of  odd  Integers  made  It  possible  to  have  an  unam¬ 
biguous  center  point,  and  It  was  this  center  point  that  was  placed  at  the  point  * 
previously  selected.  The  mapping  between  the  numbers  0  to  99  (designated  here  as 
N)  and  A  Is  given  by  Table  1.  The  distribution  of  A  given  by  Table  l  Is  baaed  upon 
the  results  ot  Rosenberg  and  Dewan. 1  While  there  remain  some  unanswered  ques- 
tlona  regarding  the  accuracy  of  this  distribution  at  this  time,  It  Is  possible  to  say 
that  tt  Is  a  reasonable  one  to  use  for  simulation  purposes.  The  actual  resulting 
value  of  KB  that  will  emerge,  however,  will  be,  Tor  present  purposes,  regarded  as 
arbitrary.  In  other  words,  the  purpose  of  this  paper  Is  to  develop  a  method  of 


lfl 


describing  turbulent  transport  In  a  stable  fluid  In  general,  and  Its  application  to 
specific  data  will  be  done  elsewhere. 


Table  1,  Distribution  of  Layer  Thicknesses  A 


N  Range 

A 

osNssa 

3Aa 

83  <  N  S  74 

BAz 

74  <  N  at  88 

?Az 

88  <  N  s  92 

BAz 

B2  <  N  S  9B 

HAz 

1  2 

In  Rosenberg  and  Dewan  and  Dewan,  some  specific  boundary  conditions 
were  used  at  the  top  and  bottom  of  the  slab.  This  caused  more  than  one  reader  to 
raise  questions,  and  It  also  gave  rise  to  confusion  about  the  mathematics  and  the 
physics  of  the  situation.  For  this  reason  we  devised  a  method  that  completely 
eliminates  any  boundary  conditions.  This  Is  the  reason  that  the  slab  In  Figure  4 
Is  located  within  an  "environment".  (The  technique  was  suggested  by  Nell 
Qroasbard,  who  was  also  responsible  for  writing  the  computer  programs  used 
below.)  The  random-number  generator  for  a  gave  values  of  uniform  probability 
from  s  ■*  o  to  a  «  max.  Layers  generated  outside  the  slab  were  not  counted  by  the 
program,  but  all  layers  within  the  slab,  Including  any  points  wtthtn,  were  counted. 
This  count  was  used  In  a  manner  to  be  described.  The  point  here  Is  that  the  layers 
outside  the  slab  served  only  to  "cause  transport"  across  the  boundaries  of  the 
latter  In  a  reasonable  manner. 

One  more  problem  remained.  The  environment  was  finite  and  therefore  It 
also  had  boundaries.  Thus,  there  was  the  queatlon  of  how  to  account  for  the 
effects  of  these  boundaries.  The  "environment"  was  chosen  to  be  much  larger 
than  the  alab,  and  the  alab  Itaelf  waa  located  In  the  center  of  the  environment,  as 
Indicated.  A  technique  to  be  described  was  used  to  trace  the  Influence  of  the  bounds 
of  the  environment.  This  Influence  propagated  toward  the  slab,  and  at  the  moment 
It  reached  the  slab  (if  It  ever  did),  the  program  would  automatically  terminate. 

The  following  procedure  was  used  to  trace  the  Influence  of  the  boundaries  at 
z  «  o  and  z  ■  max.  As  soon  as  a  layer  waa  generated  which  in  any  way  contacted 
the  boundary,  a  "cursor"  (shown  In  Figure  4  as  an  arrow)  would  be  placed  (  In 
effect)  at  the  edge  of  this  layer  which  was  nearest  to  the  slab.  When  a  subsequent 
layer  overlapped  the  cursor,  the  latter  waa  moved  again  to  that  layer  edge  nearest 
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to  the  slab.  At  the  bottom  of  Figure  4,  thla  has  occurred  twice,  and  at  the  top, 
three  times,  The  number  generator  would  generate  numbers  exclusively  for  values 

of  z  which  were  between  the  upper  and  lower  cursors  and  the  program  would  ter¬ 
minate  If  either  cursor  came  within  a  dlatance  UAz  of  the  slab.  Thua,  the  bound¬ 
aries  had  no  effect  on  the  calculation  of  flux  through  the  slab. 

In  order  to  calculate  Kg,  the  program,  In  effect,  measured  the  flux,  T,  out  oT 
the  bottom  of  the  slab,  and  used  Eq.  13),  as  was  mentioned.  aJ’/Bz  was  taken  as 
R’1,  which  was  justified  not  only  from  the  reasonableness  of  the  assumption  of 
constant  B$/8z  In  a  steady  state  process,  but  also  by  the  actual  results  to  be 
displayed  below  (see  Figure  9  D  through  1),  We  now  consider  exactly  how  7  Is 
calculated, 

As  mentioned,  F  is  the  amount  of  "material  mixed  out"  of  the  bottom  of  the 
slab  per  unit  of  time.  The  time  was  calculated  by  adding  1  to  a  register  labeled 
KOUNT  each  time  a  layer  Included  a  point  Interior  to  the  slab,  The  total  time  t 
was  calculated  from 


t  »  (KOUNT)  •  At 


where  At  Is  the  time  Interval  between  layer  formations  within  the  slab.  This 
quantity,  At,  is  to  he  regarded  as  the  average  time  between  layer  formations  In 
the  actual  situation  which  Involves  random  times  for  layer  formations,  but  the 
exact  connection  between  At  and  actual  measurements  will  be  described  below. 

The  sverage  flux  will  be  the  total  amount  of  material  (per  unit  of  uron)  to  exit  the 
bottom  of  the  slsb  (at  the  time  t)  divided  by  t.  How  Is  the  amount  of  material 
mixed  out  In  one  mixing  event  across  the  boundary  to  bo  calculated?  This  perhaps 
explained  most  simply  by  means  of  a  specific  example,  l.ot  the  thickness  be 
A  ’  5  Az,  and  let  It  be  positioned  across  the  lower  slab  boundary,  as  shown  In 
Figure  3,  l.et  the  values  of  d  «t  the  center  of  each  "cell"  of  the  layer  be  given  as 
Indicated  In  this  figure,  that  Is  by  *1  10  *8'  Only  the  two  uppermost  cells  will  be 
within  the  slab,  and  tho  rest  will  be  below,  We  wish  to  calculate  how  much  mate¬ 
rial  per  unit  of  area  Is  removed  from  the  two  upper  cells  by  the  "mixing"  or 
averaging  process. 

The  Initial  amount  of  material,  Mj,  In  the  upper  cells  (remembering  lhal  <t> 
here  will  be  a  density)  la 


M,  *  A  A*  <*4  +  *„)  (13) 

where  A  lu  the  area,  After  the  value  of  4  has  been  made  homogeneous  (within  A  by 
mixing  or,  equivalently,  by  averaging)  the  new  amount  of  mutortul  In  the  two  top 
cells,  Mf,  Is 


■-7*CT“T.'k- 
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Figure  S.  Explanation  of  the  Calculation  of  K_  by  Mean*  of  Flux  Out  of  the  Bottom 
of  the  Slab  by  Meana  of  a  "Mix-Through"  ProAaa 

Mf  ■  (AA«)(  J)  (♦1  +  *a  +  *3  +  #4  +  *8)  '  (14) 

The  amount  of  material  per  unit  area  paaatng  through  the  boundary  by  thla  mixing 
proceaa  la  given  by 


DUMP 


Mi  “  Mr 


In  general,  the  valuea  of  Mj  and  Mf  would  depend  upon  layer  thlckneaa  and  location, 
etc.  Each  time  auch  an  event  occurred,  an  amount  given  by  Eq.  (18)  waa  added  to 
a  regtater  called  DUMP.  Thus,  DUMP,  at  any  time  t,  waa  the  total  output  per  unit 
area.  The  average  flux,  F,  waa  therefore  calculated  from 


IP  ,  DUMP 

F  zrnasoRf 


..w- 


hence,  from  Eqa.  (10),  (3)  and  0$/8z  =  R'*,  we  obtain 

„  _ / DUMP •  R  \ 

Kfi-FLUX  "  ■‘  K.U11NT )  < 1 7) 

where  the  cub  Index  "FLUX"  will  distinguish  the  estimate  of  Kg  from  othere  to  be 
described  below. 

The  actual  simulation  results  from  Eq,  (l?)  will  be  shown  In  Section  0,  First, 
however,  some  other  methods  to  estimate  Kg  wtll  be  deeorlbed,  The  latter  are 
baaed  on  the  Idea  that  solutions  to  the  diffusion  equation  will  depend  on  a  constant 
value  Kg,  which  Is  to  say  that  Kg  will  be  Independent  of  Initial  or  boundary  con¬ 
ditions.  We  shall  consider  three  speolflo  cases  and  determine  If  the  estimates  of 
Kg  based  upon  continuous  theory  and  simulated  behavior  Is  self-consistent.  These 
oases  are  (a)  i(t)  Initially  la  a  delta  function  at  the  center  of  the  slab,  (b)  the 
Initial  d(s)  Interior  to  the  slab  Is  a  constant  (with  respect  to  a)  and  the  "walls"  are 
"held"  to  a  different  constant  value  higher  than  the  interior  value,  and  (c)  $(*)  la 
a  straight  line,  and  the  slab  boundaries  are  Insulated  from  the  "environment".  In 
effect,  Kg  will  be  estimated  frr  m  >U.e  diffusion  effects  In  these  three  cases, 


5.  CALCULATION  OF  KB  FROM  DIFFUSIVE  EFFECTS 

In  this  station,  we  examine  the  solutions  to  ths  continuous,  on* -dimensional, 
diffusion  equation  for  the  purpose  mentioned  above, 

8,1  Delta  Function 

The  one -dimensional  diffusion  equation  is 

K  $  ■  If  «ll» 

The  solution  in  the  case  where  the  Initial  d  distribution  la  a  delta  function  Is  well 
known  (see  Carrier  and  Pearaon,  pp.  2fl-3410),  It  la  given  by 


d<t,  a) 


-<*-*or/«Kt) 


2'/kT 


(l») 


10,  Carrier,  O.F.  and  Pearson,  C.K.  (1678)  Parttal  Differential  Equations, 
Academic  Press,  Inc,  1  "  . 


ao 


where  *Q  designates  the  initial  location  of  the  delta  function.  The  domain  of  this 
solution  Is  -*  <  7  <  ao,  and  there  are  no  boundary  conditions.  (Equation  (IB)  Is 
the  Green's  function  when  the  Initial  condition  on  #(z)  la  specified. ) 

Given  the  behavior  of  4(*>  t)  under  these  conditions,  one  can  estimate  K  from 
the  rate  of  spread  In  a  simple  way,  that  Is  by  measuring  the  1/e  halfwidth,  re,  at 
a  given  time.  The  latter  Is  obtained  from  Eq.  (IB)  by  setting  the  exponent  equal 
to  -1.  Solving  for  K 


Kb-del/a 


(20) 


(Note  that  the  "decay, "  as  used  here,  Is  spatial,  not  temporal. )  Thus,  K  can  be 
estimated  by  the  amount  of  spread  of  the  Initial  delta  function,  A  refined  estimate 
could  be  obtained  by  appropriately  combining  several  such  measurements  (taken 
at  a  number  of  times),  This,  however,  was  not  done  here.  Figure  0  shows 
schematically  the  typical  behavior  to  be  expected,  that  (s,  the  function  decreases 
In  height  and  It  spreads  as  time  Increases. 


* 
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S.2  Intnuion  from  Boundaries 

Next  we  consider  the  case  where  we  have  a  slab  such  that  4(a.  o)  ■  o  tnatde 
and  4(s,  o)  ■  40  (a  constant  larger  than  aero)  outatde.  We  aaaume  Infinite  con¬ 
ductivity  outside,  and  we  want  to  know  4(e,  t)  Inaide  (that  la,  the  development  of 
4  Inaide  aa  a  Amotion  of  time).  Notice  that,  aa  wao  the  caae  for  Figure  0,  the  alab 
la  turned  aldewaya,  with  a  going  to  the  right  In  Figure  7,  Figure  7(a)  ahowa  the 
Initial  condition  while  Figure  7(b)  ahowa  achematioally  how  4  dtffuaea  Inwards 
from  the  boundary.  The  exact  solution  of  Eq.  (18)  for  the  aeml-lnftnlte  alab  caae 
will  be  found  In  Spiegel,  p,  322.  11  (For  early  tlmea,  In  our  form  of  simulation, 
the  Intrusions  from  the  two  sides  will  not  Influence  each  other  at  all,  and  therefore 
this  solution  will  suffice  for  our  needa. )  This  solution  tat 


4(*.  t)  *  40U  ■  erf(y))  (31) 

which  In  Figure  7  would  apply  only  to  the  left  aide.  The  definitions  of  y  and  erf(y) 
are 


/  "  (3>flG) 


(23) 


and 


(23) 


A  tabulation  of  the  valuea  of  erf(y)  la  to  be  found  In  Abramowtta  and  Stvgun,  p,  311,  ** 
With  this  formulation,  one  can  start  with  a  given  4(e,  0  (from  the  simulation)  and 
estimate  K  numerically  In  the  following  manner.  Flrat,  y  la  eotlmated  via 
Eq.  (21),  from  a  given  4(n,  t)  for  a  particular  choice  of  both  a  and  t.  Then  solving 
Eq.  (22)  for  K 

kb-int  ■  ( f )  (ay)  (a4> 


11.  Spiegel,  M.  R.  (190ft)  Theory  and  Problems  of  Laplace  Transforms. 

Hcnaum  Publishing  Cc£  1 

13.  Abramowltc,  M.  and  Stegun,  I.  A.  (1964)  Handbook  of  Mathematical  Functions, 
N.B.S.  Appl.  Math.  Series  #56.  - 
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S.3  IniuUted  Slab— Initially  Linear 

Figure  8  depicts  a  slab  with  initial  distribution  of  i  given  by 

^•U  (26) 

where 

4 

a  ■  (28) 

and  the  thickness  of  the  slab  Is  given  by  R  and  maximum  4  by  4g,  We  assume  that 
the  slab  Is  Insulated  at  the  boundaries,  therefore 


ll**o,  *«r 

These  boundary  conditions  are  explained  in  any  book  on  heat  conduction,  The 
analytic  solution  of  Eq.  (18)  with  these  conditions  is  (see  for  example,  Dew  an3) 


*<«,»>■  *0  I  +  ij  l  tjjfUe  e-K<fl/R)a‘ 


where  g(n)  Is  defined  by 


0  n  even 


■  1  n  odd 


Define  the  symbol  Aj(t)  to  be  the  amplitude  of  the  first  harmonic  of  the  solution  of 
Eq.  (28).  Thus 


Aj(t)  *  (-0.40B)  e-K("/R)  * 


This  Is  obtained  from  Eq,  (28)  by  setting  z  *  o,  calculating  4/»^  •  0, 400,  regard¬ 
ing  the  factor  (1/2)  as  the  zeroth  harmonic  and  setting  b0  '  1  for  convenience  from 
this  point  on. 

Suppose  that  we  are  given  $<z,  t)  (from  the  numerical  simulation),  and  we  wish 
to  calculate  K  from  this  solution.  The  procedure  we  shall  uae  la  to  calculate  the 
flrat  harmonic  numerically,  and  to  thua  ascertain  A  at  a  given  value  of  t.  Then 
K  can  be  estimated  by  solving  Eq,  (30)  for  Kt 


f  INTERMEDIATE 
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(31) 


KB-INSUL  = 


We  now  have  the  means  to  estimate  the  bulk  diffusion  parameter,  Kg,  directly 
from  solution  behavior  for  the  three  special  cases  of  Interest.  We  now  turn  to  the 
actual  simulations  of  the  random  mixing  layer  situation, 


6.  CONNECTIONS  OF  CERTAIN  PARAMETER  VALUES  TO  EXPERIMENTAL 
DATA  AND  SIMULATIONS  OF  MIXING  LAYER  TRANSPORT 

6.1  ReUtioni  Detweon  Parameter*  At  and  A  and  Poialble  Meaiurement* 

Simulations  In  general  can  be  performed  only  for  numerically  specific  cases, 

In  this  subsection,  we  wish  to  justify  that  a  certain  choice  for  At  and  A  lit  rea¬ 
sonable  one  In  connection  with  atmospheric  studios.  We  also  wish  to  show  how 
these  parameters  could  be  measured  In  practice. 

How  could  one  measure  A,  the  turbulent  layer  thickness  ?  One  method  that  has 

in 

been  employed  (Barat  )  uses  a  hlgh-reflolutlon  wind  fluctuation  measuring  device 
mounted  on  a  balloon.  When  the  balloon  descends  through  a  turbulent  layer,  the 
turbulent  fluctuations  can  be  detected,  and  the  Information  telemetered  to  a  ground 
station.  Such  measurements  made  In  the  stratosphere  Indicate  that  a  typical  value 
of  A  would  be  around  200  m. 

An  alternative  measurement  technique  estimates  A  from  the  Richardson  num¬ 
ber  profiles  (Rosenberg  and  Dewan1),  The  present  model  was  originally  created 
for  exactly  this  purpose,  (In  other  words,  we  needed'®  method  to  convert  R{  pro¬ 
files  Into  estimates  of  transport,)  The  Rtchardson  number  profiles  are  ascertained 
from  measurements  of  the  horizontal  winds  as  a  function  of  altitude,  In  conjunction 
with  the  vertical  profile  of  the  temperature  as  well.  Triangulated  rocket-laid  smoke 
trails  Is  the  method  of  choice  for  hlgh-resolutlon  wind  profiles.  The  estimates  of  A 
can  be  obtained  from  the  fact  that  when  the  Richardson  number,  Rj,  Is  less  than  0,  25, 
the  layer  will,  presumably,  become  turbulent.  When  account  Is  taken  of  subse¬ 
quent  spreading  of  the  layer,  A  can  be  estimated.  According  to  Rosenberg  and 

Dewan,  *  A  200  m  Is  a  reasonable  number  as  estimated  by  this  technique,  For 

2 

further  discussion,  see  Dewan. 

The  parameter  At  (which  stands  for  "time  between  onsets  of  layer  formation" 
in  the  simulations  to  be  given  below)  Is  measured  less  directly,  and  Is  tied 
methodologically  to  the  method  of  estimating  A  .  In  the  case  where  A  Is  measured 

13.  Barat,  J,  (1875)  Etude  Expertmentale  de  la  Structure  du  Champ  de  Turbulence 
dans  la  MoyennTSTrato^phere,  C.ft.  Acad.  Sc.  Faria,  SBO.  Ser.  ft 


from  R,  profiles,  the  following  method  Is  appropriate.  First  (see  Rosenberg  and 
Dewan,  and  Dewan,2)  we  make  the  simplifying  assumption  that,  If  over  a  certain 
height  region  (slab)  there  are  usually  several  layers  detected,  then  such  a  profile 
will,  on  the  average,  be  replaced  by  another  profile  at  a  time  At^  later,  The  sub- 
index,  f,  la  for  "time  frame",  Tho  simplification  consists  of  replacing  the  random 
times  of  layer  formation  by  Atf  between  profiles.  Once  Rj  <  0.  23,  the  layer  to  be 
considered  la  unstable!  however,  a  certain  time  must  elapse  before  the  Instability 
Is  replaced  by  an  actual  "breakdown"  and  subsequent  turbulence.  Shortly  after 
the  latter  occurs,  Rj  will  rise  again  above  the  threshold  value  (to  a  value  of  0,  7, 
perhapa),  and  the  turbulence  will  begin  to  decay.  At  this  point,  the  R^  profile 
would  no  longer  Indicate  the  subsequent  generation  of  the  turbulent  layer.  Based 
upon  the  measurements  of  Browning1^  of  Kelvln-Helmholta  billows  In  the  atmos¬ 
phere  by  means  of  radar  and  balloon  measurements,  It  was  estimated  (Rosenberg 
and  Dewan1)  that  Atf  -  1500  s,^  In  order  to  ascertain  At  from  Atf,  one  must 
first  estimate  the  average  number  of  layers  to  be  found  at  any  time  In  the  slab, 
designated  below  as  n.  Then 


If,  for  example,  P*  designates  the  average  fraction  of  the  slab  of  thickness  R  that 
Is  turbulent,  and  suppose  we  are  given  X,  the  average  turbulent  layer  thickness, 
then 


from  P*  ■  (nX  /R). 

In  the  case  where  the  profile  of  turbulence  Is  estimated  by  means  of  In-sltu 
measurements  of  velocity  fluctuations,  Atf  would  correspond  to  the  average  time 
for  the  average  layer  to  decay  (from  time  of  measurement)  to  the  level  where  It 
would  be  below  the  "threshold  intensity"  and  would  be  called  laminar,  Based  upon 
this  given  threshold,'1  P+  would  be  ascertained  (as  well  as  A),  and  again  At  would 
be  estimated  from  Eqs.  (32)  and  (33),  For  further  discussion,  see  Dewan.2® 

^Actually,  Atf  would  correspond  to  the  average  time  between  the  observation  of 
R|  '<  0.  25  and  breakdown.  The  Utter  would  depend  on  growth  at  time  of  observa¬ 
tion  and  Rp 

■tjf  the  threshold  were  high,  P*  would  be  measured  as  small,  but  so  also  would  At; 
(because  the  layers  would  disappear  more  rapidly  below  threshold).  We  assume 
here  that  the  ratio  (P+/Atf)  remains  Independent  of  the  thresold  over  a  reasonable 
range  or  the  latter,  ^ 

®In  this  simulation  JW  «  277  m,  whereas  X  3  240  m. 

14,  Browning,  K.  A.  (1067)  Structure  of  the  atmosphere  In  the  vicinity  of  large 

amplitude  Kelvln-Helmholta  blllowB,  Roy.  Met,  Soc,  Quart,  J,  B7i2B3-2&8. 


6.2  Simulation  Resulti 


6.2.1  Kb.flttx 

Figure  9  (A<^:)  shows  the  value  of  ^.pLUX’  *?Med  on  E<**  (l7),  ■*  R  function 
of  time,  Aa  the  time  baae  (number  of  "oounta"  or  layer  mixing  eventa)  Inoreasea, 
Kg -PLUx' converIe*  to  ®  fixed  value,  Parametera  were  aaalgned  aa  followa, 

R  ■  400  polnta,  and  As  ■  50  mj  where  R  «  20  X  10®  m,  Alao,  A  ■  240  m  and 
V  X2  ■  277,  aa  determined  from  the  dlotrlbutlon  designed  Into  Table  1,  From 
Rosenberg  and  Dewan'a*  result  of  At.  ■  1500a,  we  use 


At 


ft) 


/1B00\ 

loy; 


360a 


(34) 


where  P*  ■  0,05  la  taken  from  the  same  reference,  The  simulation 
KB-FL10C  " 

Figure  0  (D— I)  ahowa  the  Important  result  that  d(s)  la  eaaentlally  unchanged 
during  thla  simulation, 


A 


_ _ _ _ . _ 
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Figure  9.  Bulk  Dtffuslvlty,  Kw  (aa  Calculated  In  the  Simulation 
Via  the  Flux)  aa  a  Function  of  Time,  (a— Time  la  given  In 
terms  "mix  through"  events.  100  B,  000  mixing  event.  Parts 
d  through  I  show  the  Initial  and  subaeqent  profiles  In  order  to 
demonstrate  that  the  mean  slope  of  d  remains  essentially  con¬ 
stant,  (The  profiles  are  shown  at  Intervals  of  5,  000  mixing 
events) 
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Figure  8.  Bulk  Dlffualvtty,  Kn  (aa  Calculated  In  the  Simulation 
Via  the  Flux)  aa  a  Function  of  Time,  (a-c)  Time  ta  given  In 
terma  "mix  through"  events.  100  «»  8,  000  mixing  event.  Parte 
d  through  l  show  the  Initial  and  aubaequent  profllea  In  order  to 
demonatrate  that  the  mean  alope  of  6  remalna  eaaentlally  con- 
atant,  (The  profllea  are  ahown  at  Intervale  of  5,  000  mixing 
event  a) 


* 


Figure  9.  Bulk  Dlffuslvlty,  K  (o«  Calculated  In  the  Simulation 
Via  the  Flux)  as  a  Function  of  Time.  (a~c)  Time  Is  given  In 
terms  "mix  through"  events.  100  ~  8,  000  mixing  event,  Parts 
d  through  I  show  the  Initial  and  subsequent  profiles  In  order  to 
demonstrate  that  the  mean  slope  of  4  remains  essentially  con¬ 
stant.  (The  proflloa  are  shown  at  Intervals  of  5,  000  mixing 
events) 
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Figure  8.  Bulk  Dlffualvlty,  K«  (as  Calculated  In  the  Simulation 
Via  the  Flux)  aa  a  Function  of Time,  (a-c)  Time  la  given  In 
terma  "mix  through"  eventa,  100  -  8, 000  mixing  event.  Parts 
d  through  l  ahow  the  Initial  and  subsequent  profiles  In  order  to 
demonstrate  that  the  mean  slope  of  8  remains  essentially  con¬ 
stant,  (The  profllea  are  shown  at  intervals  of  5,  000  mixing 


a.  a,  a  Kb.delta 

Figure  10  (A  through  H)  (compare  the  conttnuoue  ca*e)  *how*  the  re«ult*  a 
■equence  of  layer  mixing  event*  upon  an  Initial  4(a,  o),  which  la  a  delta  function, 
tn  thta  numerical  almulatlon  It  la  given  by  aero  for  all  point*  In  the  range,  R  (of 
400  point*),  except  for  the  center  point  (*Q),  which  wa*  *et  to  the  maximum  value 
(100).  At  the  arbitrary  t  KOUNT  «  30. 4X  103  (30, 400  layer  formation*),  ahd 
with  360  aecond*  aaalgned  to  each  count,  Kg.DELTA  WM  ^rom  E<1* 

The  value  of  the  peak  wa*  3. 0,  having  atarted  from  100,  and  Ita  (1/e)  value  la 
therefore  1, 1.  The  value*  of  the  end  point*  of  the  two  r#  (one  per  *ldo)  were  esti¬ 
mated  from  the  printout.  They  occurred  at  a  ■  point  1SS  and  point  273.  The  center 
wa*  located  at  point  300.  Thu*,  ualng  As  •  SO  m,  the  average  i'f  wa*  determined 
to  be  r0  ■  (3700)  meter*,  thua 

KB-DELTA  *  (4W3<flfl)°3^;iTO  *  (33) 

which  agree*  with  K^.p^uy’ 


Figure  10.  Simulated  Evolution  of  Delta  Function  When 
Expoaed  to  Random  Mixing  Layera.  A  3-D  plot,  H,  sum¬ 
marise*  the  evolution 
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Figure  10.  Simulated  Evolution  of  Delta  Function  When 
Expoaed  to  Random  Mixing  Layer*,  A  3-D  plot,  H,  turn 
marlee*  the  evolution 
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Figure  10,  Simulated  Evolution  of  Delta  Function  When  Exposed  to  Random 
Mixing  Laye re,  A  3-D  plot,  H,  summarizes  the  evolution 


8,3,3  KB.1NTR 

Figure  11  (A  through  D)  showa  the  sequence  for  the  Intrusion  case,  At 
KOUNT  *  14  X  10  ,  the  value  of  z  was  arbitrarily  chosen  to  correspond  to  point 
number  81,  or  81  X  (50  m/pt.)  *  2,  88  X  10®  m,  from  the  origin.  From  the  printout, 
4  •  0, 187  at  that  point  and  time,  The  value  of  t  was  determined  from 

t  -  KOUNT  (360)  (36) 

or  t  »  5.04  X  10°  s.  From  Eq.  (31)  (*Q  ■  1) 

crf(y)  ■  1  -  0, 167  *  0.  833  (37) 

Using  the  tables  In  Abramowtta  and  Stegun,  12  we  found  y  «  0,  98.  Thus,  Tram 
Eq.  (34) 

kb-int  ■  "-^7/  (W5$)  ■  °'M  «» 
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Figure  11.  Simulated  Evolution  of  the  Case  of  Infinitely  Conducting  Walls 
and  Diffusion  Inwards  by  Means  of  Random  Mixing  Layers.  ll-E  Is  a  sum¬ 
mary  3-D  plot 


6,2,4  KB-INSUL 

Figure  12  (A  through  Q)  shows  the  simulation  for  the  case  of  the  Insulated  slab. 
At  KOUNT  »  4  X  104  (t  ■  1,44  X  107  a),  a  Fourier  transform  was  obtained  from  the 
simulated  solution  by  means  of  a  Fast  Fourier  transform  analysis.  This  gave 
A^t)  «  (-0,381),  Inserting  this  Into  Eq.  (31) 
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Figure  IS,  Simulated  Evolution  of  "Innulated  Slab"  Case  Vie 
Random  Mixing  Layers 
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Figure  12.  Simulated  Evolution  of  "Ineulated  Slab"  Citt  Via 
Random  Mixing  Layer  a 


fl.2.6  Kg.C^LC 

We  now  uae  Eq,  (4).  Here  we  U«e 

At. 

Atb  ■  w 

(Dewan3).  Uelng  P*  »  0.05,  Atf  «  1500,  //V*'  *  277  m 


KB-CAL,C 


»•«  ”1/« 


<4 


(4 


6.8  Summary  of  Simulation  Raiulta 

Table  2  eummarlaoa  the  above  reaulta,  and  It  ehowa  In  a  convincing  manner 
that  Kg,  no  matter  which  way  It  la  eatlmated,  gtvea  a  aelf-conalatent  reuult, 


Table  2,  Distribution  Parameter,  K^,  for  Various  Cases 


KB-FLUX 

-  0.31  m2/ 

'a  -  {8.04  X  10"3)  A5/ At 

KB -DELTA 

■  0,31 

«  (8,04  x  10‘3)  A2/ At 

kb*intr 

>  0.  34 

•  (6,83  X  10‘3>  A2/ At, 

KB-1N8UL 

-  0,33 

•  (B.38X  to'3)  A2/ At, 

kb -CALC 

•  0.  33 

-  <8,  35  X  10'3)  A2/ At, 

1.45  X  10'3  (A2/ At) 
1.45  x  10*3  (A2/  At) 
1.5BX  10'3  (A2/ At) 
1.  SO  X  10’3  (A^/At) 
1.50X  10*3  (A2/ At) 


The  lMt  column  la  rendered  In  unite  of  A  and  Atf,  Thla  "nondlmenalonal"  form 
allowa  one  to  apply  the  iormaltam  to  flulda  with  other  acalea  (for  example,  the 
ocean),  Notice  that  only  single  valuea  of  a  and  t  were  uaed  to  "measure"  the 
apread  In  each  caae,  t  being  choaen  arbitrarily.  Obviously,  more  pairs  could  have 
been  uaed  to  Inarease  the  accuracy  of  Kg  estimates  from  the  diffusion  affects.  The 
consistency  ohown  In  Table  2  la  therefore  rather  remarkable,  and  we  think  one 
need  go  no  farther  at  thla  point  to  justify  the  assertion  that  the  "mixing  layer"  Kg 
of  thla  paper  la  self-consistent,  at  least  to  the  first  significant  digit. 


7.  SUMMARY  AND  CONCLUSIONS 

In  thla  paper,  a  new  diffusion  parameter  or  "effective  dtffualvlty"  has  been 

proposed,  which  la  designed  for  the  extremely  Inhomogeneous  case  of  layered 

turbulence  In  atratlfted  flulda,  Random  occurrence  In  altitude,  time,  and  layer 

thickness  were  assumed,  as  well  aa  total  mixing  within  layers,  The  last  aasump- 

2 

tlon  can  be  relaxed,  If  necessary,  aa  shown  In  Dewan, 

The  present  diffusion  parameter  describes  a  process  which,  In  essence,  Is  a 
sort  of  stochastic,  finite  difference  simulation  of  the  diffusion  equation  that  Is 
presumed  to  occur  In  the  stratified  fluids  to  be  found  In  nature.  While  the  basic 
assumptions  remain  to  be  vcrlHed,  it  la  hoped  that  the  formalism  developed  here 
will  be  useful  In  the  estimation  of  vertical  transport  by  turbulence  In  the  strato¬ 
sphere  end  upper  ocean. 


*ThU  assumes  that  P*  has  a  specific  value,  namely  P*  *  0.05. 
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Appendix 

Continuous  Ittlmate  of  Amount  of  Material  far  Unit  Aral 
PaNlng  Through  Bottom  of  "Slab”  Qua  to  an  Avaraga  Mixing 
Ivant  and  Modification  of  tha  Bitlmata  of  Kg-CALC 

Lot  the  bottom  of  the  slab  bo  located  nt  2  Zjj.  Let  2fj  designate  the  centre 
of  n  mixing  layer,  and  we  shall  assume  that  the  layer  overlaps  the  location  of  Z^. 
The  latter  condition,  where  A  Is  the  thickness  of  the  mixing  layer,  Is 

a|z0  -  ZB|  <  A  (AD 

Let  H(Zo)  be  the  height  by  which  the  mixing  layer  extends  above  Z^,  Then 

■  [(4  *  z„)  ■  *n] 

Define  M(Zo)  to  be  the  amount  of  material  per  unit  areu  which  pusses  through  the 
Z  =  Zg  boundary  as  n  result  of  mixing.  Then 

M(Z())  ■  M,(Zn)  -  Mf(Zn)  (A3) 


where  Is  defined  as  the  Initial  amount  of  mnterinl  per  unit  areu  In  the  part  of 
the  layer  whleh  Is  above  Z„,  and  Mr  Is  this  amount  after  the  mixing. 


Aaeumlng,  aB  In  the  text,  that  *<Z)  la  linear  and,  for  convenience, 
i(Z)  ■  a  Z 

where  a  la  constant,  wa  hava  (uetng  a  a  d*/Aa) 

Mi(2o)  .  (halght  of  layer  above  Zg)  X  (average  Initial  concentration) 


(A4) 


B+H  /  \ 
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dZ 
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(A  8) 


Similarly 

Mf(ZQ)  •  (height  of  layer  above  Zg)  X  (average  final  concentration) 


HIV 


*0*4 

/  ,  (&■  *)  ,z 
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(A  8) 


When  M(Z  )  le  averaged  over  all  poaalbllltlea  conelatent  with  Eq«  (Al),  we 


have 


Z„  + 


MZJ  .  i  /  (m,(Z0)  -  Mf(Z0))  dZ0 

ZB  ‘  £ 


(A7i 


After  the  calculotlona  are  carried  out, 

,3 


(AB)* 


*Thla  reault  waa  flrat  obtained  by  T.  VanZandt  (private  communication), 
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.1 


The  estimate  of  Kg  given  aa  Kq,  (4)  In  the  text  would  thus  be  changed  some 
what.  Given  that  Atb  1  tf/P*,  and  Identifying  the  flux  ?  In  Eq.  (3)  with 


-  WT) 

F  . 


(A9) 


thua,  we  obtain  the  new  eatlmatet 

(Al0) 

Thle  dlffera  from  Eq.  (4)  by  a  , 'actor  of  3/3  but,  In  aplte  of  the  agreement 
between  the  simulations  and  Eqi  14),  Eq.  (A  10)  would  be  theoretically  more 
aoaurate  than  Eq.  (4).^  Equation  (A  10)  la  presumably  the  correct  estimate  to  use 
In  future  analysts  of  experimental  data  along  these  lines)  however,  no  explanation 
Is  available  at  present  as  to  why  Eq,  (A  10)  does  not  agree  with  the  simulations, 
whereas  Eq.  (4)  of  the  text  does,  A  6*polnt  numerical,  finite-difference  formula¬ 
tion  of  this  appendix  waa  carried  out  and  the  result  was  essentially  Identical  to 
Eq.  (AS).  This  remains  to  be  clarified. 


^Equation  (4)  was  derived  by  using  Mf(Z0)  *  0  (see  Rosenberg  and  Dewan1),  This 
means  that,  Instead  of  a  "mix  through"  process,  there  was  a  "dump  through" 
process. 
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